Phase diagram of Rydberg atoms in a nonequilibrium optical lattice 



o 

(N 



C/2 

■4— > 

a- 



O 

o 



in 
in 

O 



X 




Jing Qian, Guangjiong Dong, Lu Zhou and Weiping Zhang 
Quantum Institute for Light and Atoms, Department of Physics, 
East China Normal University, Shanghai 200062, People's Republic of China 

We study the quantum nonequilibrium dynamics of ultracold three-level atoms trapped in an 
optical lattice, which are excited to their Rydberg states via a two-photon excitation with non- 
negligible spontaneous emission. Rich quantum phases including uniform phase, antiferromagnetic 
phase and oscillatory phase are identified. We map out the phase diagram and find these phases 
can be controlled by adjusting the ratio of intensity of the pump light to the control light, and 
that of two-photon detuning to the Rydberg interaction strength. When the two-photon detuning 
is blue-shifted and the latter ratio is less than 1, bistability exists among the phases. Actually, this 
ratio controls the Rydberg- blockade and antiblockade effect, thus the phase transition in this system 
can be considered as a possible approach to study both effects. 
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Introduction: Rydberg atoms with principal quan- 
tum number n 3> 1 have exaggerated atomic proper- 
ties including strong dipole-dipole interactions and long 
radiative lifetimes These properties are attractive 
for quantum information processing EL Q and quantum 
many-body dynamics simulation [J, Q . Most of these 
researches require a negligible spontaneous emission for 
reducing the quantum decoherence [fj. However, a re- 
cent research shows that when spontaneous emission is 
significant, quantum nonequilibrium dynamics could be 
demonstrated using Rydberg atom gases |7j. The sys- 
tem investigated is like a spin- 1/2 particle system, which 
undergoes a phase transition from a spatially uniform 
phase to an antiferromagnetic phase by tuning the laser 
frequency as shown in Ref. [7j. Moreover, the nonequi- 
librium induced by the spontaneous emission leads to an 
oscillatory phase. Further research has found collective 
quantum jump between a state of low Rydberg popula- 
tion and a state of high Rydberg population [8|. These 
researches open up a new window to use Rydberg atoms 
for quantum nonequilibrium dynamics simulation [9] . 

When two atoms are close, the excitation of one atom 
is prohibited by an already excited neighboring atom 
due to the level shift by the strong dipole-dipole interac- 
tion at a short distance 10]. This phenomenon is called 
dipole blockade effect and has played an imp ortant role 
in current Rydberg atom researches [lll - ll3} . However, 
recent theoretical and experimental investigations have 
shown that in a three-level two-photon Rydberg excita- 
tion scheme, there also exists antiblockade effect due to 
the Aulter-Townes splitting induced by the lower transi- 
tion [3 Ell • The possibility of enhancing the antiblock- 
ade effect by trapping atoms within a lattice has been 
proposed [lj|. So far, the exploration of physical effect 
of the coexistence of blockade effect and antiblockade ef- 
fect on quantum dynamics of Rydberg atom gases in this 
two-step excitation scheme still remains at its early stage. 

In this paper, we propose to study the quantum 
nonequilibrium dynamics of three-level atoms trapped 
within a lattice via a two-step excitation. Adopting the 
mean field approach used in Ref. Q, we predict tran- 



FIG. 1: Atoms trapped in an optical lattice are uniformly ex- 
cited to the Rydberg state |r) from the ground state \g) via 
an intermediate state |e). The lattice is divided into two sub- 
lattices, marked respectively in white and black colors. The 
neighboring atoms are coupled through the Rydberg inter- 
action V . ilp and Sl c are two Rabi frequencies respectively 
corresponding to two optical transitions. 



sition between uniform phase, antiferromagnetic phase, 
and oscillatory phase. Compared to the two-level system 
investigated in Ref. |7|, we have more control knobs, in- 
cluding the ratio of intensity of the pump light to the 
control light and that of two-photon detuning to the Ry- 
dberg interaction strength, to control the phases. When 
the latter ratio is less than 1 and the two-photon detun- 
ing is blue-shifted, bistability exists among the phases. 
The latter ratio controls the competition of the block- 
ade with antiblockade effect. Thus, both effects could be 
investigated through the phase transitions. 

Model: Our scheme is shown in Fig. [T] Ultracold 
atoms are trapped within an optical lattice. We assume 
the lattice depth is deep enough, so that the center-of- 
mass motion of the atoms can be neglected. We also 
assume there is exactly one atom per site. All the atoms 
are uniformly excited by two laser beams to a Rydberg 
state. The laser-atom coupling scheme is shown in the 
right side of Fig. [TJ where \g), \e) and |r) denote the 
ground, intermediate and Rydberg states, respectively, 
for the j-th atom, and fi p ( c ) is the Rabi frequency for the 
pump (control) laser. For convenience, we divide the lat- 
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tice into two sublattices: the black sublattice and white 
sublattice, as shown in Fig. [T] We label the black/white 
site with j/k. 

Similar configuration is used in Ref. Q, where the in- 
termediate state |e) is adiabatically eliminated for a large 
detuning of the \g) — ¥ |e) transition from pump laser fre- 
quency, and an effective two-level model can be adopted 
[16j[ ■ I n contrast, in our paper, a resonant pump laser 
is assumed to couple the ground state \g) and the inter- 
mediate state |e). Such a three- level model allows both 
the Rydberg-blockade effect and the antiblockade effect 
to be studied simultaneously. Whereas, in the effective 
two-level model, only the Rydberg-blockade was studied 

0- 

In the interaction picture, the Hamiltonian Ht de- 
scribing such a lattice of atoms can be given by (h = 1) 
Ht = Ei H j +VJ2 {jk) \r) (r\j ® \r) (r\ k . Here, H 3 de- 
scribing the interaction between the j-th three-level atom 
and two lasers is given by Hj = —5 \r) (r\- + (fl p \g) (e| ■ + 
|e) (r\j+h.c.) with 5 the two-photon detuning between 
states \g)j and \r)-. V is the dipole-dipole interaction 
between the two Rydberg atoms described by a van der 
Waals (vdW) potential. 

The quantum dynamics for this system is governed by 
the master equation of its density operator p [l7j : 



p = -i[H T ,p]+ 1 Y(\g) (e|,p|e) (g\ 



i\e) (e\ i)P } 
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(1) 

Here, decoherence due to the spontaneous emission of 
the intermediate state \e) - with its linewidth 7 is phe- 
nomenologically introduced. Dissipation from the spon- 
taneous emission of Rydberg state \r)j is ignored due to 
its relatively long lifetime (~ 100/xs) compared to 7" 1 
1- 

Now we first investigate the steady-state solutions for 
Eq. U) using the mean field approximation 18]. With 
this approximation, for the j-th atom, the Rydberg in- 
teraction \r) {r\j®J2k \ r ) ( r \k $ sums over an the nearest 
neighboring sites of the j-th site) can be replaced by a 
formula of density operator |r) (rL Efe Pk,m thus the mo- 
tional equations for the density operator elements Pj. X y 
{x,y = g,e,r) are 

p jrr = 20 c Im 

Pj,er 

(2a) 

pj.ee = 2Q C (xlYO-Pj.ge ~ lmpj <e r) ~ Pj.ee (2b) 



Pj>9 e — iQ>c (X (l Pj,rr ^Pj : ee) H~ Pj,gr^) 
Pj.er — ^^c(Pj,ee Pj,rr JCPj,gr) ^ef f Pj,< 



Pj:9 r ~ (Pj,ge XPj,er) ^effPj,gr 



Pj,ge 



(2c) 



Pi, 



(2d) 
(2e) 



In deducing Eqs.©, we have scaled the frequencies and 
time respectively with 7 and 7 _1 . We have also intro- 
duced the ratio of the pump Rabi frequency to the control 



Rabi frequency, x = ^p/^c, and the effective two photon 
detuning 5 e ft = S — V J2k Pk,rr- Without the Rydberg- 
Rydberg interaction V, Eqs. ([2]) are the typical optical 
Bloch equations. The role of the Rydberg-Rydberg inter- 
action here is to compensate the two-photon detuning S, 
as shown in the formula for 5 e ff, and thus an antiblock- 
ade effect [3, EH can be induced when 5 e ff ~ 0. 

We now start to investigate the stationary solution of 
Eqs.©. These stationary solutions can be put into three 
classes: 

(a) Uniform phase (uni): The atoms are uniformly ex- 
cited over the whole lattice, i.e., — p s k rr . Hereafter, 
the superscript s stands for the stationary solution. 

(b) Antiferromagnetic phase (AF): One sublattice has 
higher excitation than the other, i.e., Pj, rr 7^ Pk.rr- 

(c) Oscillatory phase (osc): Rydberg population os- 
cillates periodically in time between two sublattices. In 
this case, one cannot find any stable steady-state AF so- 
lutions. 

We first study the uniform phase {p S jt k \ xy — P% y )- 
From the stationary solution of Eqs. ([2]), we obtain a cubic 
equation for p s xy : 
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For AF phase, two sublattices are Rydberg-populated 
differently, and labeled 1 and 2 in the subscript of p s ., k -, . 
From Eqs.©, we obtain equations involving p\ rr and 

P2,rr : 
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(5) 



The solutions p\, 2 -, rr of the coupled equations (g]) and 
([5]) correspond to two roots of the following equation 



(p s rr ) 2 - 2Ap s „ + B = 0, 



(6) 



where A = S/V+2n 4 cX 2 (l + X 2 ) I 40* (l + x 2 ) 2 + S 2 (l + 

b = s 2 /v 2 + 4ni (1 + x 2 ) 2 / [v 2 (1 + 8 X 2 nl)] . 

When A 2 — B < 0, the quadratic equation © has 
complex roots. In this situation, the AF and oscillatory 
phases don't exist, and only uniform phase exists. In 
other cases, the quadratic equation ^ has one pair of 
real solutions. When the real solutions are stable against 
small perturbation, the system is in the AF phase; oth- 
erwise, the system is in oscillatory phase. 

Phase Diagram: In this section, we show the numer- 
ical results for the stationary solution with the experi- 
mentally tunable two-photon detuning S for x = 5 (Fig. 



3 




0.8 
0.6 

0.2 



(a) 


1 / 

If 


J 





-2 2 4 
8 (units of y) 




-2 2 4 
8 (units of y) 



FIG. 2: (color online) Stationary solutions p a rr as a function 
of detuning S with V = 5, fi c = 0.1 for different ratios (a) 
X = 5 and (b) x = 20. The stable and unstable uniform 
phases respectively show in solid line and dotted line. AF 
phase and oscillatory phase are repsectively shown in solid 
line with cross symbols, and dashed line. 



[2^,) and 20 (Fig. [2}d). Stable and unstable uniform phases 
governed by Eq. are respectively plotted in solid lines 
and dotted lines. AF phases and oscillatory phases from 
Eq. ([5]) are respectively plotted in solid lines with cross 
symbols and dashed lines. 

In Fig. HJa), there exists only one stable uniform phase 
for small detunings 5 < S c ps —0.5,. The critical point 5 C 
can be determined by A 2 = B. When 8 crosses the crit- 
ical point, the uniform phase becomes unstable, giving 
rise to a pair of AF phases. This is the pitchfork bifurca- 
tion of supercritical type [l!| , in which the uniform fixed 
point loses its stability and two stable AF fixed points 
appear. As 5 increases further to the value larger than 
1.22, the AF phase becomes unstable and the system en- 
ters into the oscillatory phase. 

The phases can be controlled experimentally by adjust- 
ing the ratio of the pump Rabi frequency to the control 
one, X- With a very strong pump transition (fl p 3> J7 C ), 
the Rydberg-excitation rate can not catch up with the 
decay rate from intermediate state, so that population 
in Rydberg state hardly changes. In Fig. [2{b) with 
X = 20 much larger than that in Fig. dja), the stable 
AF fixed points always exist and there is no oscillatory 
phase. There is a supercritical pitchfork bifurcation at 
5 = 4.67 where the uniform phase becomes stable again. 

Figures [2ta) and (b) also show the signals of Rydberg- 
blockade effect and antiblockade effect. When S = 0, 
\Pi, rr ~ P2,rr\ ~ 1 f° r AF phase. This means a reso- 
nant two-photon detuning leads to an entire Rydberg 
excitation in one sublattice, whereas the other sublat- 
tice is strongly blocked. When 6 increases to 5.0, p\ rr — 
P2 rr ~ 1-Oj an d the two-photon detuning 6 is compen- 
sated by the Rydberg energy shift (S e ff ~ 0), indicating 
the antiblockade effect. 

Since the stationary solutions are dependent on both 
the two-photon detuning 6 and the ratio of the Rabi fre- 
quencies for pump and control lasers x 5 we plot the phase 
diagram in the parameter space of (5, \) in Fig. [3] The 
number of stable uniform phases is denoted in a bracket. 
The position that <5 e // = has been marked. When 
$eff ^ (S > V J2 k Pk.rr), the Rydberg interaction plays 
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FIG. 3: (color online)Phase diagram within (S, \) space for 
V = 5, fi c = 0.1. 



negligible role, and thus the system behaves like a non- 
interacting particle system. Consequently, there is only 
one stable uniform phase, no matter how x varies. When 
S is comparable to VJ2 k pu,rr, there is complex phase 
distribution within (<5, x) space. AF phase contains a 
pair of stable nonuniform fixed points. In the regions 
marked with "uni(2)", "uni(l)+AF" and "uni(l)+osc" , 
bistability takes place between two uniform phases, or 
uniform and antiferromagnetic phase or uniform and os- 
cillatory phase. When the two-photon detuning 5 is res- 
onant (5 — 0), the system is in AF phase. When 5 is red 
shifted (S < 0) and near-resonant (S « 0), there could 
be AF phase or uniform phase depending on the ratio xi 
however, if 6 is further red-shifted, there is only uniform 
phase. 

Rydberg- Rydberg interaction: The phase diagram anal- 
ysis has shown that the Rydberg-Rydberg interaction 
strength could have strong influence on the dynamics of 
Rydberg gases. In experiments, the Rydberg-Rydberg 
interaction strength can be controlled with the atom- 
atom separation R. Now we take an isotropic vdW in- 
teraction V = —dCe/'yR 6 (d is the lattice dimension) for 
short separation R as an example to investigate the con- 
trolling of the Rydberg gas dynamics with the lattice 
size [lo,[2l|. For 87 Rb atoms, we choose 5s, bp and Ryd- 
berg state respectively as \g)j, |e) • and \r)j. The interac- 
tion coefficient for Rydberg atoms is Cq = — 870kHz/im 6 
when \r) ■ = 2is 1 / 2 [22|]j and linewidth 7 from state |e) • 
is 7/2-7T = 1.7 MHz (23|. Using these parameters, the 
normalized vdW potential is V « 1/R 6 (R is in units of 
jiim, d = 2). In this case, we obtain stationary solutions 
p s rr with S = as a function of spacing R as shown in 
Fig. @]for x = 5 (dashed lines) and 20 (solid lines). 

When R <C 1.0/im, the system is dominated by a pair 
of AF phase with one sublattice completely Rydberg- 
excited (p s rr ~ 1) and the other blockade(p^ r « 0). This 
is a demonstration of the Rydberg blockade phenomenon 
under strong Rydberg interaction [24| . From experimen- 
tal view, R <C 1.0/im can be realized in solid systems, 
but it is hard for atoms trapped within a conventional 
optical lattice [25|. However, subwavelength optical lat- 
tices can be facilitated with a recent progress in sub- 
wavelength focusing, and have been realized recently for 
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FIG. 4: (color online) Stationary solutions p" r as a function of 
lattice spacing R for \ — 5 (dashed lines) and \ — 20 (solid 
lines). 
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FIG. 5: (color online) Phase diagram in the space of initial 
conditions (pi irr , P2,rr) with (a) \ = 5, S = 2 and (b) x ~ 20, 
5 = 4. (c) and (d) show typical Rydberg population dynamics 
of oscillatory phase (p? jrr = 0.4, P2, rr = 0.6) and AF phase 
{Pi,rr = 0.6, p2,rr = 0.8), respectively. 



matter wave Bloch oscillation [26[. Thus, observing AF 
phase is promising in lattices created with subwavelength 
focusing techniques. 

When R increases, the uniform phase develops and one 
sees the two-phase (uni and AF) coexistence in the regime 
around R ~ l/.tm. After crossing over this regime, the AF 
phase disappears and the uniform phase dominates with 
R increasing further, due to the negligible value of V for 
large atom-atom separations. 

Dynamics: The dynamical evolution of the system is 
explored via a straightforward time-integration of the 
mean- field equations @. According to the phase dia- 
gram of Fig. [3l the system is in the bistable state between 
the uniform phase and oscillatory phase for \ = 5, 5 = 2. 
The results of dynamical simulation in Fig. [5fa) gives 
the information about which phase the system will finally 



evolve into by choosing different initial Rydberg popula- 
tions (pi rr , p2 rr )- If both p\ rr and p% rr are smaller than 
0.4, the final state remains in the uniform phase. By in- 
creasing p\ rr or P2. rr , due to the bistability effect, the 
system enters into oscillatory phase with Rydberg popu- 
lation periodically oscillating between two sublattices(see 
(c) for p\ rr = 0.4, p% rr = 0.6). Fig. fflb) is obtained 
in the same way except with \ = 20, 5 = 4. Under 
this set of parameters, the system is in the bistable state 
between uniform phase and AF phase according to Fig. 
[3] The uniform phase now becomes dominated. This 
is because with a far-off resonant two-photon detuning 
and a weak control laser, the transition from intermedi- 
ate state to Rydberg state becomes very difficult, leading 
to small Rydberg population and in turn weak Rydberg 
interaction. Therefore the system tends to stay in the 
uniform phase. However, as long as p\ rr and p\ rr are 
both significantly large, the system will still evolve into 
the AF phase. Fig. [S^d) presents a typical Rydberg pop- 
ulation dynamics in the case of AF phase for p\ rr = 0.6, 
p\ rr = 0.8. 

Conclusions: We have studied mean- field phase dia- 
grams in a lattice system of three-level atoms with two- 
photon Rydberg excitation and spontaneous emission- 
induced dissipation. We have shown that the system has 
rich phases, including antiferromagnetic phase, uniform 
phase, and oscillatory phase. We can control the phases 
using either the ratio of the intensity of the pump to the 
control light, or that of two-photon detuning to the Ryd- 
berg interaction strength. The latter ratio actually tunes 
the competition of the Rydberg-blockadc effect and the 
antiblockade effect. We found when the Rydberg interac- 
tion strength V is much larger than two-photon detuning 
S (the Rydberg-blockade dominates), the system can be 
in the AF phase. When V is comparable to 6, the com- 
petition of the Rydberg-blockade eand the antiblockade 
effects leads to bistability of the system. When \S\ is far 
larger than V, the system is in the uniform phase. Our 
work shows the possibility of studying antiblockade effect 
within a lattice (l5j . 
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